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BEAMFORMED  NEARFIELD  IMAGING  OF  A  SIMULATED 
CORONARY  ARTERY  CONTAINING  A  STENOSIS 


1.  INTRODUCTION 

Coronary  artery  disease  (CAD)  is  a  primary  precursor  of  heart  attacks,  the  leading 
cause  of  death  in  the  United  States.  CAD  is  the  buildup  of  plaque  in  the  coronary  arteries 
resulting  in  a  condition  referred  to  as  stenosis  in  which  blood  flow  is  restricted  and  the  oxygen 
supply  to  the  heart  muscle  is  decreased.  The  end  result  of  arterial  stenosis  is  frequently 
myocardial  infarction  and  subsequent  necrosis  (i.e.,  permanent  heart  muscle  damage). 
Presymptomatic  detection  of  coronary  stenosis  before  heart  damage  has  occurred  is  extremely 
important  because  the  disease  could  be  treated  with  diet  and  drugs  rather  than  with  coronary 
bypass  surgery  or  angioplasty.  Early  CAD  diagnosis  on  a  wide  scale  would  save  the  U.S. 
health  care  system  billions  of  dollars  and  would  allow  a  much  better  quality  of  life  to  those  who 
suffer  from  heart  disease. 

The  most  common  method  to  measure  blood  flow  and  stenosis  in  coronary  arteries  is 
by  an  invasive  test  called  an  angiogram.  During  this  procedure,  a  catheter  injects  a  radiopaque 
substance  into  the  coronary  arteries,  while  an  x-ray  screen  simultaneously  reveals  the  arterial 
blockages  of  the  patient  to  the  doctor.  However,  an  angiogram  is  costly  and  places  the  patient 
at  significant  risk.  For  many  years,  a  nonintrusive  method  called  ultrasound  has  been  used  by 
cardiologists  to  image  the  coronary  valves.  This  technique,  which  exploits  the  compressional 
wave  field  in  the  human  body,  typically  uses  an  actuator  to  emit  low-level  energy  into  the  body 
where  it  is  absorbed  and  reflected  to  varying  degrees.  An  accompanying  array  of  sensors 
receives  the  energy  back  and  an  image  is  formed  with  a  volume-scanning  focused  beamformer. 
Although  ultrasound  can  also  image  the  inside  of  a  near-surface  artery,  such  as  the  carotid 
artery,  it  cannot  image  a  coronary  artery  due  to  its  depth,  its  smaller  size,  and  the  nearby 
scattering  of  energy. 

In  this  report,  a  new  method  that  is  proposed  to  detect  and  localize  arterial  stenoses  uses 
the  shear  wave  space-time  field  of  a  solid  body  and  phased  array  technology  to  image  the 
underlying  volume.1  The  approach  is  based  on  the  turbulent  normal  wall  forces  that  are  present 
in  and  downstream  from  a  stenosis.23  These  wall  forces,  while  extremely  small,  exert  a 
pressure  on  the  wall  of  the  artery,  which  causes  a  displacement  in  the  surrounding  soft  tissue. 
This  displacement  occurs  in  the  form  of  two  elastic  waves  in  the  medium:  a  compressional 
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wave  and  a  shear  wave.  At  low  frequencies  (100  to  500  Hz),  the  shear  wave  dominates  the 
response  and  can  be  measured  at  the  surface  of  the  solid  by  an  array  of  sensors.  Then, 
nearfield  (focused)  beamforming  imaging  algorithms  allow  detection  and  localization  of  the 
stenosis.  An  in  vitro  experiment  is  included  to  validate  this  method  and  corresponding 
simulations  are  presented. 
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2.  FOCUSED  BEAMFORMING 


The  beamformer  is  defined  by  assuming  that  x(CO0)  is  a  zero-mean,  stationary, 
complex  vector  of  N  random  sensor  outputs  that  are  Fourier  transformed  at  fixed  frequency 
(Oq.  The  image  scalar  pixel  value  of  the  conventional  focused  beamformer  (CFB)  at  Cartesian 
coordinates  x-v  yv  and  z-  is  expressed  by 


c{xbybzba0) 


d(x^  ’  ’  6>0 )  R(6>q  )  d(x;' ,  ,  Zp  tOq ) 

d(*i , y(- ,z;-,£00)H  d(x;- , yf , z,- , £t)0 ) 


(1) 


In  equation  (1),  R(o>0)  =  E[x(fi)0)  x(tf)0)H]  is  the  iV-by -N  cross-spectral  density  matrix 
(CSDM),  E[  ]  is  the  statistical  expectation  operator,  and  d(jq,y,-,Zj,tf>o)  is  die  N- by-1  focused 
imaging  vector  with  the  nth  row  expressed  as 


dn(xi  »  ■)/  >  "-i  )  ft  ®XP 

te(«)r 

where  /?  is  the  geometric  spreading  loss  factor,  c  is  the  average  wave  speed,  and  rj-(n)  is  the 
apparent  Euclidean  distance  from  the  ith  image  point  to  the  nth  sensor  located  at  x(n),  y(n),  and 
z(n).  This  Euclidean  distance  is  expressed  as 


~j  COpn(n) 


(2) 


n(n)  =  *J[xi  -  x(n)]2  +  [yf  -  y(n)]2  +  \zt  -  z(n)]2  . 


(3) 


In  addition,  the  image  scalar  pixel  value  of  the  reduced  variance  distortionless  response 
(RVDR)  focused  beamformer  at  Cartesian  coordinates  x-,  y  •,  and  z;-  and  at  fixed  frequency  £»0 
is  defined  by 


Mxi,yi,Zi,co0) 


w(x;,y;-,zt-,6)Q)H  R(ft)p)  w(x;-,y„z;-,a>o) 

W(Xj ,  yi ,  Zi  >  Oq  )  ^  Zj  >  ©o ) 


(4) 


where  w(x,-,y,-,z,-,fi)o)  *s  die  iV-by-1  focused  minimum  variance  distortionless  response 
(MVDR)  sensor  array  filter  vector.  This  vector  is  expressed  as 


w(xf , yf ,Z;,<00)  =  R(Cd0 )  1  d(x;- , yt , z; , <o0 )  , 


(5) 
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A 

where  d(xIsy,-,Zi,©o)  *s  given  in  equation  (2)  and  R(fi)0)  is  the  diagonal  enhanced  CSDM 
defined  by 


R(dO0)  =  R(®o) +  e 


Trace  [R(ty0)]] 

1  *  J 


(6) 


where  e  is  the  scalar  enhancement  factor  and  I  is  the  TV-by-TV  identity  matrix.  It  is  noted  that  if 
e  equals  infinity,  equation  (4)  becomes  a  conventional  beamformer,  and  if  e  equals  zero, 
equation  (4)  becomes  an  MVDR  beamformer.4 
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3.  SIMULATION 


Numerical  simulation  of  this  problem  consists  of  analytically  formulating  the  CSDM, 
inserting  it  into  equations  (1)  and  (4),  and  then  calculating  the  simulated  beamformed  image. 
The  analytical  model  of  the  turbulent  flow  energy,  based  on  previous  laboratory  wall  pressure 
measurements  downstream  of  a  blockage,3  is  modified  to  include  energy  in  the  region  of  the 
blockage.  Additionally,  a  spatial  correlation  term  is  added  to  the  model  to  account  for  the 
spatial  coherence  of  the  turbulent  flow.5  The  expression  for  the  streamwise  wall  pressure 
autospectrum  is  given  as 


P(x,co0 )  =  p(x)ex  p 


ry\x\(o0^ 

U. 


exp 


U, 


(7) 


where  p{x)  is  the  amplitude  of  the  wall  pressure  shown  in  figure  1,  x  is  the  streamwise  location 
of  the  flow  measured  from  a  reference  point  (m),  y  is  the  spatial  decay  constant 
(dimensionless),  and  U c  is  the  fluid  convection  speed  (m/s).  In  figure  1,  the  ordinate  is  the 
root  mean  square  (rms)  wall  pressure  nondimensionalized  by  the  fluid  density  p,  the  square  of 
the  constriction  velocity  Up  the  constricted  diameter  d,  and  the  unconstricted  diameter  D;  the 
abscissa  is  the  spatial  location  x.  For  this  study,  the  spatial  decay  constant  was  -0. 125  and  the 
fluid  convection  speed  was  the  flow  velocity  through  the  constricted  portion  of  the  tube. 
Equation  (7)  is  used  to  form  the  M-by-M  CSDM  E(&)0)  of  the  source  (blockage  and 
downstream  energy),  which  has  the  element  in  the  Mth  row  and  vth  column  expressed  as 


Euv(o) o)  = 


PiXuMXy) 


{max*  [/>(*)]} 
p(xu)p(xv) 


J  y\xu-xv\o)0)  f  j  co0(xu -xv)^ 
2  CAr  rr  CAr 


{max*  [p(*)]}' 


'eXP 


Uc 

r|*«-*v  K 

Ur 


Ur 


exp 


ijOQ  Ov~*M) 

Ur 


U<V 


V  <  U 


(8) 


For  the  source  model  used  here,  the  indices  u  and  v  range  from  1  to  21  and  correspond  to  the 
circled  data  points  in  figure  1,  and  maxx[p(x)]  is  the  maximum  value  of  the  pressure  with 
respect  to  spatial  location  x.  Once  this  matrix  is  known,  the  simulated  CSDM  can  be 
determined  using 


R(o0)  =  &s 


N 


Trace  [d(q>0)  E(©0)  Dh(co0)]J 


D(m0)  E(o)0)  DH(<y0)  +  GqI  ,  (9) 
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where  D(®0)  is  the  N-by-M  matrix  from  the  nth  sensor  to  the  nth  source  whose  nnth  element  is 
expressed  as 


o)  = 


U(n,u)f 


exp 


— j<a>o  l{n, 


(10) 


where  l(n,u )  is  the  Euclidean  distance  from  the  nth  sensor  to  the  nth  source,  given  as 


l(n,u)  =  ^[x(rf)  -  x(n)]2  +  [y(n)  -  y(n)]2  +  [z(n)  -  z(n)]2  ;  (11) 

Oq  is  the  spectral  density  of  the  zero-mean  spatially  uncorrelated  additive  noise  at  each  sensor; 
and  <r2  is  the  signal  at  the  peak  of  the  spatial  spectrum. 
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4.  EXPERIMENT 


A  homogeneous  urethane  mold  containing  a  latex  surgical  tube  was  manufactured  for 
the  in  vitro  tests,  as  shown  in  figure  2.  The  surgical  tube,  intended  to  simulate  the  mid  left 
anterior  descending  coronary  artery,  has  an  inner  diameter  of  0.00318  m  (0.125  in.)  and  is 
embedded  at  a  depth  of  0.0318  m  (1.25  in.)  in  the  urethane.  The  urethane,  which  is  Hexcel 
195-RE,  represents  the  surrounding  soft  tissue  and  has  a  shear  wave  speed  of  7  m/s  at  100  Hz 
and  12  m/s  at  1000  Hz.  The  100-  to  1000-Hz  shear  wave  speed  is  approximated  by  using  a 
linear  dispersion  relationship  with  respect  to  frequency.  (It  should  be  noted  here  that  the 
slowest  shear  wave  speed  available  in  a  commercial  grade  urethane  is  somewhat  faster  than  soft 
human  tissue  shear  wave  speed  at  low  frequencies.6)  The  surgical  tube  contains  a  flow 
blockage  (nylon),  which  represents  a  coronary  (partial)  occlusion  having  a  length  of  0.0127  m 
(0.5  in.)  and  a  hole  with  an  inner  diameter  of  0.00159  m  (0.0625  in.),  as  shown  in  figure  3. 
This  blockage  was  chosen  because  it  has  a  50-percent  diameter  reduction  (or  a  75-percent  area 
reduction),  which  is  similar  to  the  reduction  used  by  other  researchers  who  have  measured 
break  frequencies,7  velocity,8'10  and  pressure31112  caused  by  a  constriction  in  a  cylindrical 
tube.  The  in-flow  end  of  the  nylon  occlusion  contains  a  45-degree  taper  to  provide  a  relatively 
smooth  entrance  region. 

Water  is  gravity  fed  into  the  surgical  tube  from  a  holding  tank  and  captured  below  the 
urethane  mold  in  a  calibrated  flask,  as  shown  in  figure  4.  Although  this  is  a  steady-state  flow 
experiment,  it  has  been  shown  that  the  force  levels  are  nearly  identical  to  a  corresponding 
pulsatile  flow  experiment.12  This  technique  permits  measurement  of  the  flow  velocity  from 
which  the  corresponding  Reynolds  numbers  based  on  obstructed  and  unobstructed  diameters 
can  be  computed.  The  moving  fluid  transmits  a  relatively  small  normal  force  into  the 
embedded  tube,  which,  in  turn,  transmits  energy  into  the  urethane.  This  energy  propagates  to 
the  surface  of  the  mold,  primarily  in  the  form  of  shear  waves,  where  it  produces  a  local 
displacement  field  that  fluctuates  in  space  and  time.  A  polyvinylidene  fluoride  (PVDF)  sensor 
array  is  used  to  measure  this  field  and  the  beamformer  is  used  to  localize  the  energy  emissions 
in  the  three-dimensional  solid.  The  experiment  was  conducted  in  an  electromagnetically  (EM) 
shielded  anechoic  chamber  to  minimize  ambient  noise  contamination. 

All  of  the  measurements  in  this  report  were  taken  with  a  two-dimensional  array  that 
consists  of  15  sensors  (elements)  on  two  orthogonal  axes,  as  shown  in  figure  5.  A  common 
middle  sensor  is  shared  by  the  nine  sensors  (equally  spaced)  on  the  primary  axis  and  the  seven 
sensors  (not  equally  spaced)  on  the  secondary  axis.  The  element-to-element  midpoint  spacing 
on  all  of  the  primary  axis  sensors  is  0.00953  m  (0.375  in.).  The  element-to-element  midpoint 
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spacing  of  each  of  the  outer  three  sensors  on  the  secondary  axis  is  also  0.00953  m  (0.375  in.), 
and  the  distance  from  the  center  sensor  to  its  adjacent  sensor  is  0.0286  m  (1.125  in.).  Each 
group  of  three  outboard  sensors  on  the  secondary  orthogonal  axis  is  built  into  a  separate  box 
that  is  attached  to  the  main  box  by  hinges  so  that  the  array  can  conform  to  nonplanar  surfaces. 
A  coordinate  system  is  chosen  such  that  the  primary  axis  is  the  x-axis  and  the  secondary  axis  is 
the  z-axis.  The  y-axis  is  the  image  depth  of  the  urethane  block.  The  coordinate  system  origin 
is  the  intersection  of  the  primary  and  the  secondary  array  axes. 

The  sensing  elements  of  the  array  are  PVDF  film  strips  that  are  0.0381  m  (1.50  in.) 
long  and  approximately  0.00927  m  (0.365  in.)  wide.  Attached  to  the  middle  of  each  strip  is  a 
nearly  circular  piece  of  plastic,  which  has  a  0.00927-m  (0.365-in.)  approximate  diameter  and  a 
0.00159-m  (0.0625-in.)  thickness.  This  plastic  contact  piece  is  the  sensing  area  of  the  element 
and  is  used  to  facilitate  an  omnidirectional  response  at  the  short  wavelengths  that  occur  at  the 
slow  wave  speeds  and  corresponding  frequencies  involved  in  this  study. 

The  array  elements  were  connected  to  a  20-channel  low-noise  preamplifier  with 
lightweight  multistrand  coaxial  cables  (Cooner  Corporation).  The  signals  were  then  passed  to 
a  bank  of  amplifiers  (Ithaco  Model  451)  with  the  low-frequency  rolloff  set  to  1000  Hz  and  then 
to  a  multichannel  programmable  filter  (Precision  Filter)  with  the  high-frequency  cutoff  set  to 
800  Hz.  Finally,  the  signals  were  passed  to  a  data  acquisition  system  (LabView)  residing  in  a 
computer  (Apple  Power  Macintosh  8100/100)  where  the  time  series  data  were  stored  for  post 
processing. 
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Figure  4.  Laboratory  Configuration 


10 


Figure  5.  Two-Dimensional  Array  With  1/8-Inch  Surgical  Tube 
Depicted  as  a  Dashed  Line 
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5.  RESULTS 


Two  different  geometrical  configurations  were  numerically  modeled  and  corresponding 
experimental  measurements  were  conducted  to  validate  the  models.  The  first  is  an  ideal  case 
where  the  main  axis  of  the  array  is  parallel  to  the  surgical  tube  and  the  middle  sensor  is  directly 
over  the  occlusion,  as  shown  in  figure  6.  Figures  7,  8,  9,  and  10  are  comparisons  of  the 
simulation  and  the  data  for  200,  300, 400,  and  500  Hz  with  the  CFB.  Figures  1 1, 12,  13,  and 
14  are  comparisons  of  the  simuladon  and  the  data  for  200,  300,  400,  and  500  Hz  with  the 
RVDR  focused  beamformer.  All  the  surfaces  are  normalized  such  that  the  maximum  value  is  0 
dB.  The  experimental  data  were  computed  with  48  fast  Fourier  transforms  that  were  each  1024 
data  points  long;  the  sampling  frequency  was  2016  Hz  per  channel,  which  corresponds  to  a 
frequency  bin  resolution  of  1.969  Hz.  The  wave  speed  used  to  compute  all  the  images  was 
linearly  interpolated  for  each  specific  frequency  with  a  value  of  7  m/s  at  100  Hz  and  12  m/s  at 
1000  Hz.  This  resulted  in  a  wave  speed  of  7.56  m/s  at  200  Hz,  8.1 1  m/s  at  300  Hz,  8.67  m/s 
at  400  Hz,  and  9.22  m/s  at  500  Hz.  The  spreading  parameter  /3  was  1,  which  represents 
spherical  spreading  at  each  source  location  and  the  RVDR  enhancement  factor  was  2.  The 
Reynolds  number  in  the  surgical  tube  (ReD)  was  1790  and  the  fluid  velocity  was  0.57  m/s. 
The  Reynolds  number  through  the  constriction  (Red)  was  3580  and  the  fluid  velocity  was  2.28 
m/s.  It  is  noted  that  the  effective  constricted  diameter  (the  vena  contracta)  is  slightly  smaller 
than  the  actual  constricted  diameter.13  Numerical  simulations  with  smaller  constricted 
diameters  produced  the  same  results  as  those  using  the  actual  constricted  diameter. 

In  figures  7  though  14,  the  x-axis  represents  the  coordinate  that  is  parallel  to  the  main 
array  axis  with  the  middle  sensor  centered  at  x  =  0.  The  y-axis  is  image  depth  into  the  solid. 
The  two  parallel  black  lines  centered  at  y  =  3.18  cm  on  figures  7  through  14  represent  the 
surgical  tube  in  the  urethane.  The  solid  black  rectangle  centered  at  x  =  0  and  y  =  3. 18  cm  is  the 
occlusion.  The  fluid  is  flowing  from  the  right  side  of  the  image  to  the  left  side.  The  yellow 
rectangles  at  y  =  0  depict  the  sensor  location  on  the  main  axis  of  the  array.  The  sensors  on  the 
secondary  axis  of  the  array  are  not  shown.  The  color  bar  on  the  right-hand  side  of  each  image 
is  the  scale  of  the  image  in  decibels.  The  minimum  surface  value  plotted  is  -10  dB. 

The  second  geometrical  configuration  is  the  case  where  the  main  axis  of  the  array  is 
orthogonal  to  the  surgical  tube  and  the  blockage  is  between  two  of  the  sensors  of  the  array. 
Shown  in  figure  15,  this  is  representative  of  a  human  test  where  an  artery  with  an  occlusion  is 
at  some  unknown  location.  Figures  16,  17,  18,  and  19  are  comparisons  of  the  simulation  and 
the  data  for  200,  300,  400,  and  500  Hz  with  the  CFB,  and  figures  20,  21,  22,  and  23  are 
comparisons  of  the  simulation  and  the  data  for  200,  300,  400,  and  500  Hz  with  the  RVDR 
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focused  beamformer.  The  Reynolds  number  in  the  surgical  tube  (ReD)  was  1820  and  the  fluid 
velocity  was  0.58  m/s.  The  Reynolds  number  through  the  constriction  (Red)  was  3640  and  the 
fluid  velocity  was  2.32  m/s.  The  simulated  and  experimental  constants  used  in  the  first  case 
were  also  used  in  this  one. 

In  figures  15  though  23,  the  x-axis  represents  the  coordinate  that  is  parallel  to  the  main 
array  axis  with  the  middle  sensor  centered  at  x  =  0.  The  y-axis  is  image  depth  into  the  solid. 
The  arrow-marked  black  circle  centered  at  y  =  3. 18  cm  and  x  =  2.38  cm  on  figures  16  through 
23  represents  the  surgical  tube  in  the  urethane.  The  fluid  is  flowing  up  and  out  of  the  image 
and  is  orthogonal  to  the  image  plane.  The  yellow  rectangles  at  y  =  0  depict  the  sensor  location 
on  the  main  axis  of  the  array.  The  sensors  on  the  secondary  axis  of  the  array  are  not  shown. 
The  color  bar  on  the  right-hand  side  of  the  images  is  the  scale  of  the  image  in  decibels.  The 
minimum  surface  value  plotted  is  -10  dB. 
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Figure  6.  Array  and  Blockage  Location  for  First  Experiment 
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Figure  7.  Simulation  and  Experimental  Data  of  CFB  for  Experiment  One  at  200  Hz 
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Figure  9.  Simulation  and  Experimental  Data  of  CFB  for  Experiment  One  at  400  Hz 
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Figure  10.  Simulation  and  Experimental  Data  of  CFB  for  Experiment  One  at  500  Hz 
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Figure  11.  Simulation  and  Experimental  Data  of  RVDR  Focused  Beamformer 

for  Experiment  One  at  200  Hz 
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Figure  12.  Simulation  and  Experimental  Data  of  R  VDR  Focused  Beamformer 

for  Experiment  One  at  300  Hz 
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Figure  13.  Simulation  and  Experimental  Data  of  RVDR  Focused  Beamformer 

for  Experiment  One  at  400  Hz 


29/(30  blank) 


Figure  14.  Simulation  and  Experimental  Data  of  RVDR  Focused  Beamformer 

for  Experiment  One  at  500  Hz 
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Figure  16.  Simulation  and  Experimental  Data  of  CFB  for  Experiment  Two  at  200  Hz 
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Figure  1 7.  Simulation  and  Experimental  Data  of  CFB  for  Experiment  Two  at  300  Hz 
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Figure  18.  Simulation  and  Experimental  Data  of  CFB  for  Experiment  Two  at  400  Hz 
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Figure  19.  Simulation  and  Experimental  Data  of  CFB  for  Experiment  Two  at  500  Hz 
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Figure  20.  Simulation  and  Experimental  Data  of  RVDR  Focused  Beamformer 

for  Experiment  Two  at  200  Hz 
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Figure  21.  Simulation  and  Experimental  Data  ofRVDR  Focused  Beamformer 

for  Experiment  Two  at  300  Hz 
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Figure  22.  Simulation  and  Experimental  Data  ofRVDR  Focused  Beamformer 

for  Experiment  Two  at  400  Hz 
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Figure  23.  Simulation  and  Experimental  Data  ofRVDR  Focused  Beamformer 

for  Experiment  Two  at  500  Hz 
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The  CFB  and  the  RVDR  focused  beamformers  each  produced  images  that  accurately 
detected  and  localized  the  blockage.  The  CEB-measured  images  were  slightly  better  than  the 
RVDR-measured  images,  especially  in  the  second  experiment.  This  is  typically  the  result  of  a 
model  parameter  mismatch  between  beamformer  parameters  and  actual  measurement  values. 
High-resolution  beamformers  (such  as  the  RVDR  beamformer)  tend  to  be  substantially  more 
sensitive  to  parametric  mismatch  than  conventional  beamformers.  Specifically,  errors  in  shear 
wave  speed,  sensor  position,  spatial  coherence  of  the  noise,  spatial  coherence  of  the  source, 
and  geometric  spreading  loss  factor  will  produce  model  mismatch  that  distorts  the  beamformed 
image  to  varying  degrees. 

The  images  have  better  resolution  as  the  frequency  increases.  This  is  a  result  of  the 
shorter  wavelength  energy  producing  a  beam  with  a  smaller  area.  In  addition,  beams  that  are 
produced  with  this  type  of  imaging  algorithm  provide  less  resolution  as  the  distance  of  the  focal 
point  into  the  imaged  volume  increases.  The  image  depth  of  3.18  cm  (1.25  in.)  was  chosen 
because  this  is  the  shortest  distance  from  the  chest  wall  to  the  left  anterior  descending  coronary 
artery  in  the  average  human.14 

It  is  noted  that  the  simulations  and  experimental  data  are  in  good,  although  not  exact, 
agreement.  There  are  several  possible  reasons  for  this  variance: 

•  The  turbulent  energy  is  modeled  with  a  line  segment  source  that  has  no  physical 
extent  in  the  radial  direction,  whereas  the  experiment  uses  a  tube  that  is  1/8  inch  in 
diameter. 

•  The  point-to-point  correlation  function  of  the  source  model  was  derived  from  a  flat 
plate  model,  whereas  the  point-to-point  correlation  of  the  tube  may  be  slighdy 
different. 

•  The  line  segment  source  model  used  was  an  rms  pressure  frequency-independent 
measurement  rather  than  a  frequency-dependent  measurement  that  could  have  a 
different  spatial  shape  at  each  frequency. 

Additionally,  the  blockage  used  in  the  previous  source  measurements3  had  a  blunt  inlet  region 
compared  with  the  tapered  inlet  region  used  in  this  experiment.  Nevertheless,  the  sidelobe 
energy  in  most  of  the  measurement  images  is  in  approximately  the  same  locations  as  the 
corresponding  numerical  simulations,  which  is  a  good  indication  that  the  source  model  in  the 
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simulations  is  similar  to  the  actual  turbulent  energy  in  the  experiment.  Finally,  the  background 
noise  assumed  in  the  numerical  simulations  is  spatially  uncorrelated  with  zero  mean  and 
uniform  variance,  whereas  the  noise  in  the  actual  experiment  could  have  angular  extent.  This 
type  of  noise  mismatch  in  the  model  would  cause  a  bias  error  for  energy  location  because  the 
focusing  process  has  inaccurately  described  the  effective  noise  spatial  extent. 
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6.  CONCLUSIONS  AND  RECOMMENDATIONS 


In  summary,  the  turbulent  flow  energy  in  and  downstream  from  the  occlusion  creates  a 
normal  wall  pressure  on  the  tube  that  is  converted  to  (predominantly)  shear  waves  in  the  solid 
at  frequencies  of  100  to  500  Hz.  These  waves  propagate  to  the  surface  where  they  can  be 
measured  by  an  array  signal  processing  system  that  includes  nearfield  beamforming  algorithms 
to  perform  imaging  of  the  spatial  shear  wave  field.  Imaging  the  volumetric  region  under  the 
array  produces  an  accurate  representation  of  the  internal  energy  field.  The  use  of  a  -3  dB 
contour  image  deflection  criterion  in  the  study  allows  the  occlusion  to  be  detected  and  localized 
to  within  1  cm  of  depth  and  crossdepth. 

Further  research  should  investigate  the  capability  of  this  method  to  more  accurately 
locate  the  blockage  and  more  fully  characterize  the  temporal  spectrum  of  energy  that  it 
produces.  Moreover,  it  is  recommended  that  embedded  stiffeners  representing  the  rib  cage  be 
added  to  the  urethane  mold.  This  will  change  the  energy  sensed  by  the  array  because  the 
waves  will  scatter  off  and  reradiate  from  high-density  objects  in  the  solid.  Such  volumetric 
inhomogeneity  will  be  far  more  representative  of  the  in  vivo  conditions  that  the  experiment 
must  replicate  to  achieve  the  ultimate  goal  of  this  effort,  which  is  the  detection  of  early  stage 
arterial  disease  in  the  human  body. 
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